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Abstract 



The Lehmann-Newton-Wu mass matrix, which was recently applied to neutri- 
nos, is further investigated. The analytic results presented earlier are confirmed 
numerically for the solar density profile of the Standard Solar Model. The combined 
analysis of atmospheric and solar neutrino data favors the MSW solution over the 
vacuum-oscillation solution. The total rates from the solar neutrino detectors favor 
two regions in the (7711,772-2,772.3) mass space, with the heaviest one, 7773 ~ 0.05 eV. 
The spectrum distortion reported by the Super-Kamiokande collaboration favors one 
of these two regions, which has two light (7771,7772 0.003 eV) neutrinos. 



1 Introduction 



Neutrino flavor transitions (e.g. v e — > v^) are the most widely accepted explanation of 
the atmospheric neutrino anomaly and the solar neutrino problem (SNP) pQ. In the most 
popular scenarios for such flavor conversions, vacuum oscillations and matter enhanced 
transitions (MSW) [2], the neutrino conversion and survival probabilities depend on two 
types of neutrino-specific parameters. One of them is the mixing matrix U defined from 
v a = ^UaiVi, where a = e,/i, or r and the sum runs over the mass states, Ui, z/ 2 , v%. 
The other is the differences of masses squared, Am^ = m\ — m 2 , where the latter refer to 
the neutrino mass states V{. There is no widely accepted mass matrix, therefore the con- 
ventional way of analyzing the results from neutrino oscillation experiments is to consider 
the mixing matrix as independent of the mass-squared differences. But when we assume a 
certain mass matrix M in the Lagrangian (in a current basis) 

£ = vMv + h.c, (1.1) 

then the mixing will explicitly depend on the masses through the diagonalization of M. 
The maximal mixing arising from the democratic texture jH] is a widely discussed example 
of this kind (see also |lj). 

The mass matrix proposed by Lehmann, Newton and Wu (LNW) for quarks [Sj, was 
recently applied to neutrinos in [6], and gross features of the model were explored for solar 
and atmospheric neutrinos. We shall here examine the implications of that model in more 
detail, focusing on solar neutrinos, and compare the vacuum-oscillation scenario with the 
MSW scenario. Also, we will for the different cases discuss the energy spectrum of the 
recoil electrons induced by solar neutrinos. Constraints from atmospheric neutrino data 
will be taken into account. 

2 The model 

Let us briefly review the model being considered [13 E]. One assumes that there are three 
neutrinos with a real mass matrix of the form: 

/ d \ 

M = d c b , (2.1) 
\ b a J 

where b 2 = 8c 2 . This symmetric mass matrix is diagonalized by the matrix U, 

M d = U^MU, (2.2) 

where Md is a diagonal matrix whose entries are the mass eigenvalues. The diagonalization 
involves solving a cubic equation for a. 

For an arbitrary set of mass eigenvalues, (mi, m2, m^), the model yields two possible 
mixing matrices (two physical solutions of the cubic equation), referred to as Solution 1 
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and Solution 2 j^j. The corresponding mixing matrices U are displayed in Fig. ^ for 
m 2 < m 3 (the physical region of the model extends a little beyond the triangle given by 
mi < itt-2 < rnz). As illustrated in those figures, the mixing elements are determined by 
the two mass ratios mi/rn^ and 1712/1713. 



Solution 1 Solution 2 




Figure 1: Absolute values of mixing matrix elements, U e i, U e 2, U^, . . . vs. mi and 7712 
(m 3 = 1) for Solutions 1 and 2. Dotted contours: 0.1, 0.2, 0.3, 0.4; dashed: 0.5, 0.6, 0.7; 
solid: 0.8, 0.9 (regions above 0.5 are shaded). Phases are chosen such that the angles in 
Eq. ()2.3|) are positive. Elements (2, 2) and (3, 1) are negative. 

The neutrino oscillation data give quite strong restrictions on both the mixing and 
the mass-squared differences, when treated as free parameters. Also, a model for the 
neutrino masses should not violate the experimental upper bounds on the neutrino masses 
[Zj: m Ue ^ 10 eV, m Vll < 0.17 MeV, m„ T < 18.2 MeV. Furthermore, it should satisfy 
the mass constraints from cosmology |Hj : m Vf ^ 94 eV, be capable of explaining the 
apparent oscillation of atmospheric neutrinos, with Am 2 ~ 10~ 3 eV 2 , and give a plausible 
explanation of the SNP, with Am 2 ~ 10 -10 eV 2 (vacuum oscillation) or Am 2 ~ 10~ 5 eV 2 
(MSW). The natural scale implied by these values for the mass-squared differences is far 
too small to be in conflict with any of the above-mentioned mass bounds. 

We denote the mass states in such a way that nil < m 3 . Of the two independent Am^, 
one will be reserved for an explanation of the atmospheric neutrino anomaly, and the other 
will be assigned to a solution of the solar neutrino problem. Thus we are trying in one 
model to satisfy the two most serious neutrino problems. According to our notation, the 
smallest mass-squared difference will be either Am^ or Am| 2 . 

In the present model one can fit directly in terms of masses as opposed to mass-squared 
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differences and mixing factors. When the values of the mixing elements are given by the 
masses, we do not expect as good a theoretical fit to the observed data as when these 
parameters are considered as independent. For given values of U and Am^ the probability 
expression for flavor conversion will be the same whether the mixing matrix is dependent 
on the masses or not. Therefore we know to some extent from the traditional analyses 
where in the parameter space we can expect acceptable fits for the model we are analyzing. 
For example, for solar neutrinos it is known that both vacuum- and matter-influenced 
neutrino transitions require small in order to give good fits. From Fig. ^ we see that 
this is achieved when both 7711/777.3 and m^jvn^ are small. And from the same figure it is 
apparent that in this region the mixing elements U e \ and U e 2 are very sensitive to small 
changes in the masses. This means that by small changes in the mass ratios, the mixing 
elements can be adjusted close to the best-fit values obtained when U and the masses are 
considered as independent. 



Detector 


Observed Rate 


Predicted Rate 


Obs./Pred. 


Homestake [5] 


2.56 ±0.16 ±0.16 


7 7+1.2 
' • ' -1.0 


0.332 ±0.029 


SAGE dm 


fi7 9+7.2 +3.5 
u ' - z -7.0 -3.0 


129l| 


n co 1+0.062 
u ' Jzi -0.059 


GALLEX [HI 


77.5 ± 6.2i|;? 


12911 


601+ 058 
u.ou±_ 0060 


Kamiokande ^2] 


2.80 ±0.19 ±0.33 


K 1 5+1.0 


0.544 ± 0.074 


Super-K [T5] 


2.44 ± 0.0512$ 


C 1 c;+1.0 
d - id -0.7 


n 474+0.020 

u -^' ^-0.017 



Table 1: Measured and predicted (no oscillations) (Hj event rates for solar neutrino detec- 
tors. The flux units for the three upper detectors are given in SNU, 1 event per second per 
10 36 target atoms. For the scattering detectors the fluxes are rated in units of 10 6 cm~ 2 s _1 . 
The rightmost column shows the fraction between observed and predicted fluxes. 



In the analysis of the solar-neutrino data (summarized in Table 1) we will see that, for 
small values of m\ and 777,2, the two mixing solutions give approximately the same results. 

In the traditional analyses of the data from neutrino-oscillation experiments one often 
expresses the amount of mixing in terms of angles. We adopt the following, common 
parameterization of the mixing matrix 



U 



C\2 C13 S12 C13 S13 

-Cl2 S13 S 23 - S12 C 23 -S12 S13 S 2 3 + C12 C 23 C13 S 2 3 
-C12 S13 C 23 + S12 S23 -S12S13C23 - C12S23 C13C23 



(2.3) 



where Ci 2 = cos6 , 12 , s 12 = sin 6*12, etc. This will be particularly useful when we discuss the 
MSW mechanism. If #13 = we get the two-generation case. 

For vacuum oscillations as well as for the MSW scenario, we shall discuss four cases, 
small Am^, small Am| 2 , together with Solutions 1 and 2. 
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3 Atmospheric-neutrino constraints 

We shall not here perform a detailed x 2 fit to the atmospheric-neutrino data. Instead, we 
will use best-fit values for U and AmL from the literature as constraints and compare those 
with the acceptable regions we find for the solar neutrino data. The experimental data 
indicate that the "missing" muon neutrinos can be explained by a close to maximal «-> v T 
oscillation which requires small U e3 . Also, the result from the CHOOZ experiment 
[TH] , interpreted in terms of three generations and our choice for the Am^- values, indicates 
small U e s. For masses where our model gives good fits to the solar neutrino data, we actually 
have C/ e 3 <C 1. This criterion allows us to neglect matter effects for atmospheric neutrinos 



With three neutrino families the oscillation probability between arbitrary flavor states 
can be written as 



where R is the distance traveled since the neutrinos were created, and E is the neutrino 
energy. 

For the case of atmospheric neutrinos, the factor R/E is quite small as compared with 
the case of solar neutrinos. Thus, the term in Eq. (|3.1|) with a small Am^ (henceforth 
referred to as the SNP value, ~ 10~ 5 eV 2 or ~ 10~ 10 eV 2 ) can be neglected. When this 
small SNP term involves Am^ we get 



where we use the approximation Am^ — Am§ 2 . From we know that the best fit 
between data and the oscillation hypothesis occurs at Am^ ~ 3.2 x 10" 3 eV 2 . This means 
that the lowest possible value of the heaviest mass state is m 3 = 0(0.05 eV). 

In Fig. El we have plotted the amplitude in Eq. ()3.2|) for Solutions 1 and 2. One should 
keep in mind that the approximation leading to Eq. ()3.2j) is only satisfied close to the 
diagonal, m\ = m 2 . For Solution 1, an unrealistically large m 2 /m 3 would be needed to get 
an acceptable — v T mixing. For Solution 2, we observe that large — v T mixing could 
be reached for values of Am^ which are marginally compatible with the solar-neutrino 
constraints. We also note that the Am^ which is required for the MSW effect will allow 
a larger amplitude for atmospheric-neutrino transitions than a mass-squared difference 
relevant for the vacuum oscillation solution to the SNP, since the MSW effect requires a 
larger mi- 

When Am 2 2 assumes the SNP value, i.e., m<i ~ 777,3, then the atmospheric — * v T 
transition probability can be written as P v ^ Vt [R) — ^U^U^ sin 2 (Am^R/AE) instead of 
()3.2)1 . where we have used the approximation Am^ ~ Am^. The maximal amplitude for 
this oscillation is about 0.5 for both mixing solutions. This value is too low to be of real 
interest. 



D2I- 




j<k 



(3.1) 




(3.2) 
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Sol. 1 , small Am 2 , 



Sol. 2, small Am 2 , 




Figure 2: Amplitude 4U^ 3 U^ 3 as function of mass ratios mx/ms and rri2/m 3 , relevant for 
atmospheric neutrinos under the assumption of small Am^. The maximal values are 
around 0.9 and 1 for Solutions 1 and 2, respectively. Contours corresponding to 0.7 and 
0.9 are indicated. 

In summary, the atmospheric-neutrino data favor parameters corresponding to Solu- 
tion 2 with Am^ <C ^ m 3i — Am 2 2 . They also favor the MSW solution for the solar- 
neutrino problem, since this mechanism requires a larger 1712, and thus a larger amplitude 
for atmospheric-neutrino oscillations. 

4 Vacuum oscillations 

Vacuum oscillations provide a plausible solution to the solar- neutrino problem ^Hl EH] ■ In 
this scenario, the length of the relevant oscillation cycle (see Eq. (jHHJ) is about one Sun- 
Earth distance, i.e. ~ 1.5 x 10 11 m. The corresponding Am^ is several orders of magnitude 
smaller than the one which is relevant for atmospheric neutrinos. As mentioned, it will 
be either Am^ or Am 2 2 which is relevant for solar neutrinos. Thus, the allowed values of 
Am§ 2 will be of the order relevant for either atmospheric or solar neutrinos. 

One necessary condition for acceptable fit between data and the vacuum oscillation 
hypothesis is that the smallest Am| ~ 10~ 10 eV 2 . The other independent mass-squared 
difference is held fixed at 3.2 x 10~ 3 eV 2 , in order to be able to accommodate the atmo- 
spheric neutrino data. To find the optimal mixing matrix we have to search through a 
large number of values for the mass ratios m\jm 3 and m 2 /m 3 , which we do by varying m\ 
and Ami . 
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4.1 Small Am^ 

The interesting region is here close to the diagonal in the m\-mi plane (see Fig.GJ. The best 
fit is achieved at mi = 3.1 x 1CT 6 eV, Am^ = 8.0 x 10 -11 eV 2 which gives \ 2 — 4.1. This 
X 2 -value is about 0.1% higher than the value obtained in an analysis where U and Am^ 
are treated as independent. This good agreement is explained by the high sensitivity of the 
mixing matrix to the mass values. Fig.|3]shows how the relevant mixing elements depend on 
mi when Am^ and Am.32 are held fixed at 8.0 x 10~ n eV 2 and 3.2 x 10 -3 eV 2 , respectively. 
As can be seen from this figure, the mixing solutions are equal for mi < 10~ 3 eV. For 
Solution 1, Ufj,s <C 1, which means no «-> v T for atmospheric neutrinos. For Solution 2 
this amplitude is ~ 0.4. This value is unchanged over a broad range of mass values, 
including the best fit value for the solar neutrino solution. In Fig. 0] we have plotted the \ 2 
values corresponding to Fig. El When mi is reduced beyond the best-fit value, the mixing 
matrix approaches the identity matrix, hence the fit becomes poor. 



0.8 
0.6 
0.4 
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10" 
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Figure 3: Dependences of U e i, U e 2 and U e z on m\. The mass-squared differences Am^ and 
Am 2 2 are held fixed at 8.0 x 10~ n eV 2 and 3.2 x 10~ 3 eV 2 , respectively. The two solutions 
are equal up to mi ~ 10~ 3 eV. For Solution 1, the mixing element U e ^ is close to zero for 
all values of mi. 



4.2 Small Am] 2 (m 2 /m 3 ~ 1) 

In this case we can not get an equally good fit for solar neutrinos as we had for small 
Am|i. It turns out that the best fit is achieved when we, in the expression ()3.1|) . have 
vanishing amplitude for the oscillation factors containing the largest AmL JHj- This can 



7 



Vacuum 



100 



Sol. 2 



80 



Sol. 1,2 




CM 




Sol. 1 



L 
10 



-6 



10 



-4 



10 



-2 



m. [eV] 



Figure 4: Fit to the solar-neutrino data for small Am^. The values for that parameter 
and for Am§ 2 are the same as in Fig. El The best fit gives x 2 = 4.1 at mi = 3.1 x 10~ 6 eV. 

not be achieved within the model, neither for Solution 1 nor for Solution 2, see Fig. ^ For 
Solution 1 the best fit gives x 2 ~ 6.1, while for Solution 2 one gets U e2 — U e3 ~ l/V%, 
with the best fit result x 2 — 5.0. Also, the relevant amplitude for atmospheric neutrino 
oscillations would be rather low in this case. 

4.3 Summary on vacuum oscillations 

If we consider only the SNP, none of the best fits can be excluded. The best-fit point is 
achieved with small Am^ and at practically the same mass values for Solutions 1 and 
2. Including atmospheric neutrino data, the preferred solution would be small Am^ and 
Solution 2. But the maximum possible oscillation amplitude for z/ M <-> u T is only 0.4, so 
the vacuum oscillation solution to the SNP is disfavored by the model in jB] . On the other 
hand, the recoil spectrum has an enhancement at high energies, as suggested by the data 
(see Fig. E|) . 



For an exponential solar density, the three-flavor MSW equations have recently been solved 
analytically for arbitrary masses [201 ■ Here, since we assume two widely separated values of 
AmL, we use the more conventional approach [21]. This allows for a numerical comparison 
of the two approaches. The agreement was found to be very good. For example, for (mi, 



5 The MSW interpretation 
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m 2, m 3) = (0.002, 0.01, 0.057) eV and E = 1 MeV, the two approaches give v e survival 
probabilities of 0.668 and 0.660. For some of the mass points we also check the exponential 
solar density approximation against the solar density profile calculated in a standard solar 
model fUE2], by performing numerical integrations on the latter. The agreement is within 
a few percent. 

In the two-generation MSW case where the angle and mass-squared difference are varied 
independently, one arrives at two possible solutions in the parameter space spanned by 
sin 2 (28) and Am 2 . There is one region with low \ 2 a ^ small angles, and another minimum 
at large angles [23J. As we will see, this is also the case for the model in p. 

For solar neutrinos, Solutions 1 and 2 are nearly identical for the masses relevant at 
the best-fit point. Therefore we will divide the presentation of our results into two parts, 
small Ami and small AmL 



5.1 Small Am^ 

For small values of m^jva^ and very small values of mi/m 3 we get [7 el ~ 1 for both 
mixing solutions. These are practically equivalent for solar neutrinos, both of them have 
a best fit with \ 2 — 1-8 at m i — 2-8 x 1CT 6 eV, Am^ = 7.0 x 1CT 6 eV 2 , where we have 
imposed Am^ = 3.2 x 10~ 3 eV 2 . The latter quantity is held fixed at this value, which 
is compatible with the atmospheric neutrino data [E]. Because it is large compared to 
the SNP-relevant Am| 1; its exact value is, within the uncertainty, unimportant for the 
solar- neutrino analysis. These mass values correspond to small mixing angle, see Eq. (|2.3j) . 
With the above-mentioned values of the mass-squared differences held fixed, we show in 
Fig. El (left panel) how U e \ and x 2 vary when mi is slightly perturbed around its best-fit 
value. 

As previously noted, for atmospheric neutrinos we consider only vacuum oscillations, 
due to the smallness of In the parameter region shown in Fig. El (left panel), U e % < 
10~ 3 . From Eq. (j3.2j) it follows that the v T oscillation amplitude is 0.05 and 0.6 for 

Solutions 1 and 2, respectively. Thus, Solution 1 is unacceptable in this parameter region. 
We also note that the \ 2 value is quite sensitive to the mass values mi and m 2 . 

In Fig. 03 (right panel) we show how the model compares with the data in the region 
known as the "large-mixing-angle solution" . In this region, which requires somewhat larger 
values of mi, the relevant mixing elements are very sensitive to changes in the mass ratios 
mi/m3 and m2/m 3 , which are still rather small. For the masses shown in this figure, U e ^ ^ 
10~ 2 . The best fit (varying m\ and m?) has \ 2 = 3-8 and occurs at vii\ = 2.1 x 10~ 3 eV, 
AmJi = 3.2 x 10~ 5 eV 2 , with the constraint Am|i = 3.2 x 10~ 3 eV 2 . For Solutions 1 and 
2 the atmospheric neutrino-oscillation amplitudes are 0.07 and 0.7, respectively. Thus, in 
the LMA region, Solution 1 is unacceptable^] 

These two low-% 2 regions can be combined into a single plot, Fig- El The solid contours 
enclose regions where x 2 < 7.8. Because we have two degrees of freedom (four different 

* Updated results from the Gallium detectors and Super-Kamiokande were recently presented at the 
XIX International Conference on Neutrino Physics & Astrophysics. Those values shift the SMA solution 
to mi = 2.2 x 10~ 6 eV with Am^i = 1-2 X 10~ 5 eV 2 . The LMA solution is relatively less changed. 
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Figure 5: Mixing element U e i and x 2 vs - m ii f° r Solution 1 (dashed line) and Solution 2 
(solid line) in the small-mixing-angle (left panel) and large-mixing-angle (right panel) MSW 
case. In the SMA case, the best fit has \ 2 = 1-8 and occurs at mi = 2.8 x 10~ 6 eV, 
Am^ = 7.0 x 10" 6 eV 2 . In the LMA case, the best fit has \ 2 = 3.8 at m x = 2.1 x 10" 3 eV, 
Am^ = 3.2 x 10 -5 eV 2 . In both cases, we have imposed Am| x = 3.2 x 10~ 3 eV 2 . 



kinds of detectors and two parameters) and the global minimum has Xmm = 1-8; from the 
relation x 2 = Xmin + 5-99, these contours corresponds to 95% C.L. In the same plot we also 
show the contours for atmospheric «-> v T amplitudes equal to 0.7, 0.8 and 0.9. No point 
in this figure gives a mass-squared difference Am^ outside of the 68% C.L. region given 

by nig. 



5.2 Small AmJ 2 (m 2 /m 3 ~ 1) 

In this case the mixing elements are far less sensitive to small changes in the mass ratios. 
And for this reason the mixing elements are less adjustable within the relevant values of 
the mass-squared differences. For Solution 1 we get a best fit \ 2 ~ 12. 

For Solution 2 it is unavoidable to get U e2 = U eS = 1/ y/2, see Fig. [TJ The predicted fit 
to the observed values corresponds to \ 2 — 15, which is quite poor. Small Am 2 2 is also 
unfavored in view of the atmospheric neutrino data. 



6 Spectrum distortion 

The Super-Kamiokande detector measures the recoil energy of the scattered electrons com- 
ing from the process v + e~ —> v + e~ [24]. For each of the solar sources, the form of the 
neutrino energy spectrum is measured in laboratories, and is independent of the total flux. 
Thus, the Super-Kamiokande measurements of the upper part of the 8 B energy spectrum 
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Sol. 2, m 3 =0.057 eV 



0.9 / / L 

LMA ->//' i 

0.8 /'/ 

0.7 ' 

^Z2S<- SMA 

10" 6 10" 5 10" 4 10" 3 
m 1 [eV] 

Figure 6: The lower left and right-hand contours correspond to the 95% CL. for the 
small-mixing-angle and large-mixing-angle solutions of Fig. respectively. The global 
minimum ("x") has Xmm = 1-8- The dashed, dash-dotted, and dotted lines indicate values 
of4^ 3 f/ r 2 3 . 

gives restrictions on the neutrino parameters, independent of solar models. If an energy 
dependent suppression of the flux should be established by such a method, it would be 
a very strong hint for neutrino oscillations because it would in principle not involve the 
uncertainties associated with solar models. However, the hep neutrino flux is poorly known 
|25| . This fact tends to reduce the significance of these measurements. 

In Fig.0we have plotted the ratio of the measured counting rate against the expectation 
in the absence of oscillations [21] . If there is no neutrino flavor transition, this ratio should 
be a straight horizontal line. However, it seems to increase for E > 13 MeV. In the 
same figure we have tried to reproduce the observed spectrum within three scenarios: the 
vacuum, SMA, and LMA solution. These are actually all acceptable, but, as concerns the 
spectrum, the vacuum and SMA solutions are favored. 



> 

CD 

CO 
I 

O 



14 
12 
10 

8 
6 
4 

o 



7 Conclusion 

We have extended earlier work on the application of the Lehmann-Newton-Wu mass ma- 
trix jS] to neutrinos by investigating also the vacuum oscillation case, and by studying in 
more detail the MSW solutions. Two new fits to the solar neutrino counting rates have 
been found, beyond those of jS], one at rri\ ~ 3 x 10~ 6 eV, with m 2 ~ 3 x 10~ 3 eV and 
772,3 ~ 0.05 eV, corresponding to a small mixing angle solution in the MSW interpretation, 
and a vacuum-oscillation solution with m\ and 777,2 of the order 10~ 6 -10~ 5 eV. Both of these 
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Figure 7: Energy- dependence of recoil electrons, due to oscillations with masses (mi, 772,2, 
777,3) as indicated in the figure: dashed and dash-dotted for LMA-MSW, solid for SMA- 
MSW, and dotted for vacuum oscillations. The data points show the observed attenuation 
w.r.t. the Standard Solar Model (no oscillations) I14j. 



solutions lead to electron energy recoil spectra that are in better agreement with the data 
than the LMA solutions, but they are only in marginal agreement with the atmospheric 
neutrino data. 

For both the vacuum and MSW solutions of the solar neutrino problem, only hierar- 
chical mass values gives viable fits. In both cases, best fits are obtained when the largest 
mass value is 777,3 — 0.05. When atmospheric neutrino constraints are taken into account, 
only the large-angle part of the MSW effect gives an acceptable fit to the observed data. 
The result from the combined analyses of the atmospheric and solar neutrino data will be 
further constrained with the forthcoming long-baseline experiments and SNO data |2f)j . 
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